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Abstract 

We investigate the phase structure of the two-dimensional lattice Gross-Neveu 
model formulated with the Wilson fermion action to leading order of 1/A^ 
expansion. Structural change of the parity-broken phase under the influence of 
finite temperature and chemical potential is studied. The connection between 
the lattice phase structure and the chiral phase transition of the continuum 
theory is clarified. 
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I. INTRODUCTION 



In spite of significant effort expended over the years , thermodynamic studies of lattice 
QCD with the Wilson fermion action have been lagging behind those with the Kogut- 
Susskind action. The origin of difficulty is the explicit breaking of chiral symmetry due to the 
Wilson term, with ensuing complications in analysis of chiral properties. It has become clear 
only recently that the finite-temperature phase diagram has an unconventional structure , 
having a region of spontaneously broken parity and flavor symmetry in addition to the 
usual parity- flavor symmetric phase. While this development has considerably clarified a 
number of puzzling features observed in numerical simulations in the past there are still 
questions needing further elucidation. One of the questions is how the continuum limit is 
to be taken with the new phase diagram. The parity-flavor broken phase has an extension 
which depends on the temporal lattice size, and the tuning of parameters necessary to 
achieve a continuum limit having chiral symmetry has not been explored in detail. Another 
interesting question is how the phase diagram generalizes for finite quark chemical potential 
corresponding to finite baryon density. As is well known, finite-density studies of lattice 
QCD have been plagued with serious problems 0. A conceptual understanding of the 
phase diagram is a prerequisite in numerical studies of this difficult problem. 

In this article we carry out an analytical exploration of these problems in the context 
of the two-dimensional Gross-Neveu model on the lattice formulated with the Wilson 
fermion action The model shares with QCD the feature of asymptotic freedom and 

spontaneously broken chiral symmetry. Furthermore it is analytically solvable in ex- 
pansion. These points make the model an useful arena for exploration of theoretical issues 
with lattice QCD thermodynamics with the Wilson quark action. Indeed the model has 
provided significant hints for understanding the structure of the finite-temperature phase 
diagram 0. 

In the continuum the phase diagram of the Gross-Neveu model on the plane of tem- 
perature T and fermion chemical potential n was determined some time ago |Tl|. In the 
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leading order of expansion, the (T, /i) plane is divided into two phases, a chirally broken 
phase at low temperatures and small chemical potential and a symmetric phase at high tem- 
peratures and large chemical potential, separated by a phase boundary. Along the phase 
boundary, the transition is of second order for small chemical potential, which however, 
changes into a first-order transition for large chemical potential. Our aim in this article will 
be, first, to determine the phase diagram of the lattice model for finite temporal lattice sizes 
corresponding to finite temperature and finite chemical potential, and, second, to study how 
the continuum phase diagram is recovered as one takes the limit of continuum space-time. 

This paper is organized as follows. In Sec. ||, after a brief review of the Gross-Neveu 
model in the continuum, we formulate the lattice model with the Wilson fermion action at 
finite temperature and chemical potential, and analyze the behavior of the effective potential 
toward the continuum limit. In Sec. |T| the phase structure of the model at zero temperature 
and zero chemical potential is studied. In Sec. |^ effects of finite temporal lattice size {i.e., 
of finite temperature) on the phase diagram is examined, and the continuum extrapolation 
is studied. The case of finite chemical potential is treated in Sec. |V] where we consider in 
detail how the difference in the order of phase transition observed in the continuum theory 
arises in the context of the lattice model. We conclude with a summary in Sec. |V1[ 



The Gross-Neveu model in two-dimensional Euclidean continuum space-time is defined 
by the Lagrangian density. 



where is an A^-component spinor field and denotes the coupling constant. Our conven- 
tion for the two-dimensional 7 matrices is 



II. ANALYTICAL EXAMINATIONS 



A. Continuum Theory 




(2.1) 



7i = '^2, 72 = 0-1, 75 = 0-3 = ^7i72- 
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For massless fermion m = 0, the model possesses U{1) chiral symmetry defined by 

^ ^ e^^'^^V' , ^ e-^^^^V'. (2.3) 

In terms of tlie bosonic fields introduced by 

9^ - 9^ - 

^ = ~]v^^ ' ^ ^ "^Y^^^sV" , (2.4) 

chiral transformation ( p.3|) represents a rotation on the a-n plane by an angle 26. 

Statistical properties of the system at finite temperature T can be examined by restricting 
the imaginary time extent of the space-time to /5 = 1/T, and replacing energy integrals for 
fermions by Matsubara mode sums over half-integer values according to 

fca ^ cj„ = ^(n + 1/2) ,neZ. (2.5) 

In order to describe the system at finite chemical potential /i, we add nip'^tp to C where i/j^i/j 
is the fermion number operator. These replacements, together with the introduction of the 
effective fields (^.41), leads to the action given by 

S = 



Vil^iOfi -m- a - «757r + fx'j2)i' - 7r^[(^ + vr ) . (2.6) 

"^9 



— oo 



2tx Jo 2tx 

To leading order of expansion, the ground state of the model is determined by the 
minimum of the effective potential 1^((t, tt) for a constant a and vr given by 

l/(a,7r) = — (a^ + vr^)-- / ^ln[(a + m)2 + vr2 + P + (a;„ + z/i)2], (2.7) 

where M represents an ultraviolet cutoff and nmax = i3M/{2tt). In the chiral limit m — > 0, 
y is a function of = + vr^. At zero temperature and chemical potential, the minimum 
of the effective potential is located at s = A = 2Me^^^^'^ ^ 0, which signals spontaneous 
breakdown of chiral symmetry. As the cutoff M is removed M — > oo, (7^ must converge to 
zero to keep A finite, that is, the model is asymptotically free. 

At finite temperature and chemical potential, Wolff |]TT| analyzed the effective potential 
and determined the phase diagram in the T-/i plane, which is reproduced in Fig. |l|. 



Inside the boundary ABC chiral symmetry is spontaneously broken, while the region outside 
of this line is chirally symmetric. The point 5 is a tricritical point separating a second-order 
phase boundary AB from a first-order one BC. 



B. Lattice Theory with Wilson Fermion Action 

Consider a two-dimensional lattice of a lattice spacing a. The lattice Gross-Neveu model 
with the Wilson fermion action is defined by [^] 

Aat = -^^[^{x) {r --iu)ij{x + jj) +i'{x) (r + 7^) ^/'(x - /i) ] 

+ ^(2r + 5ma)ij{x)ij{x) - ^[ + gl{ij{x)i-i^ij{x)f]. (2.8) 

where 5m is a mass counterterm. We take the coupling constant and g"^ for the two 
four- fermion interaction terms to be different |T^, whose reason will become clear below. 
Hereafter we let Wilson parameter r be unity. 

Similar to the continuum theory, we define a pair of bosonic fields according to 

= -^-ip-ip + , It = -j^tpi-f^ip , (2.9) 

In order to consider the system at finite temperature T, we take a lattice with = 
1 / (Ta) sites in the temporal direction. A finite chemical potential fi is introduced through 
a modification of the hopping factor in the temporal direction The Lagrangian of the 
model then takes the form, 

Aat ^ ^(^) (1 - 7i) + i) + i>ix) (1 + 7i) - i) ] 

-^IH^) (1 - 72) e-''" ij{x + 2) + ^(x) (1 + 72) e ^(x - 2) ] 

+ -(2 + 5ma)^{x)ij{x) - ^{a - 5mf - ^7r\ (2.10) 
a 2gi 

With these modifications, the effective potential Vl = Va^ in lattice units is given by 

1 1 
VLicTL, til) = ^{aL - Sm^f + ^vri - J ^ ^J{0 , (2.11) 



with 



^ Nt-I 

J{i) = — ln[A- Scos(a;„a + i//L)], uj^a 

i V'- 



T n=0 



27r , 1, 



A = 2 + {aL + 2f - 2{aL + 2) cos{ + tt^, 
S = 2((7l + 2) -2cos^, 



(2.12) 
(2.13) 
(2.14) 



where (7^ = era, ttl = 7ra, ^m^, = 5ma and ^ = /cia are quantities in lattice units. 
The Matsubara mode sum in J'{^) can be carried out as follows: 



-I Nt id 



1 



g IXl-^/Nt ^2-Kni/NT _ _^iiL+iT^ /Nt ^-^t^ui/Nt 



n=l 



= In 

Nt 



B^NT 



^ent ]J(i_e- 

n=l 



-E-iiL-i'n/NT ^'Kni/NT\ 



X 



_ ^-E+hl+itv/Nt ^-2Trni/NT^ 



n=l 



= In + + ^ In 1 1 + e^^-'^^)^^ | + ^ In 1 1 + e-(^+'^^)^^ | , 
2 Wj" iVj' 



where 



E = i 



cosh-\A/B) iiA/B>0 
cosh-^\A/B\+i7r HA/B < , 



and we have used the formula 



N 



n(i+ 



xe 



2nni/N\ 



1 - i-x) 



N 



n=l 



(2.15) 



(2.16) 



(2.17) 



The real part of E has the meaning of the energy of the fermion one-particle state for a 
given value of , t^l and dimensionless space momentum ^. 



C. Continuum Limit 

The lattice Gross-Neveu model as we defined above explicitly breaks chiral symmetry 
due to the Wilson term. Toward the continuum limit a — > 0, an examination of effects of 
symmetry breaking is possible through an expansion of the effective potential in powers of a 
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p!0| . We briefly recapitulate the analysis here as it raises an important point for the study 
of phase diagram in the following sections. 

We consider the effective potential at T = 0{Nt = oo) and /i = 0, and write it as 



VLiaL, ttl) = ^{aL - Stul)' + - £^ ^ In [Do(0 + 0,(0] , (2.18) 



where 



DoiO = E sin'^. + E (1 - cose.) + ^1 + (2.19) 

Di(0 = 2aL E (l-cos^-)- (2-20) 

The continuum limit of the effective potential is determined by terms of 0(0^) with ai and 

tcl regarded as 0{a). Since Di{^) = 0{a) while -Do(0 = ^(1); make an expansion in 

terms of -Di(0 ( |2.18| ) obtaining 

1 1 d?^ 

VL(aL, ttl) = — ((Tl - (5mL)2 + — vri - / — — ^lnDo(0, 

-^9 a ^9tx •'-'^ l^^J 

+ V — / — ^ (2.21) 

h ri (27r)2Do(0" ^ ^ 



The integral over lnDo(0 ( 2.21|) can be estimated by an expansion around o"|+7r| = 0, 



paying attention to logarithmic divergence of the first derivative of the integral with respect 
to a\ + tt\. The result reads 

, /.N ^ / 2 , 2^ 1 / 2 , 2x, + ^T? 



where 

Co = 0.220634 ■■• . (2.23) 



In the sum over n in (|2.21| ) only the terms up to and including n = 2 contribute to O{o?), 



and a\ + 7r| in -Do(0 can be set to zero in this limit. We then find that 



The integrals ( |2.24| - |2l25| ) clearly show that the Wilson term distorts the effective potential 



in the a direction relative to that of tt. Collecting the results together we find 



+3^-('^i + 'ri)l"^i^ + 0(a»), (2.26) 
47r e 



The necessity for introduction of the two couplings and should now be clear [jlO 
chiral symmetry would not be recovered toward the continuum limit unless one tunes the 
couplings. The limit of massless fermion further requires a tuning of the mass parameter 
5mL to remove the linear term in a. A natural tuning will be 

^ = 4 + 4C2 + 0(a), (2.27) 

and 

^-^ = -2C, + 0{a'). (2.28) 
Let us introduce the A parameter through 



A = -e-"/^' , c = e2^^°-^^^2 _ 0.57160 ■ ■ • . (2.29) 
a 



With the tuning (|2.27|) this corresponds to a coupling renormalization given by 

^ = C„-2C, + ^luJ^. (2.30) 

with which the effective potential (in physical units) takes the standard continuum form 

1 2 2 

V{a, n) = -{a' + rr') H^^) + 0{a). (2.32) 

It is straightforward to extend the above analysis to the case of finite temperature and 
chemical potential. No new divergences arise in these cases, and the same tuning of pa- 
rameters gl,gl and 6mL as specified in (p.27H2.28|) is necessary for the correct continuum 
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limit. Hereafter, we often use this tuning ignoring the 0{a) corrections. This is equivalent 
to taking 

T{gl) ^ {-2C,gl , {l/gl + ^C,)-') (2.33) 

as a tuned point on the {drriL^g'^) plane. 

We learn from the present analysis that understanding of the continuum limit requires 
an elucidation of the phase diagram of the lattice model in the three-dimensional parameter 
space of (gf^, g1, Stul). This will be our basic viewpoint in the following sections. 

III. PHASE DIAGRAM AT T = AND /x = 

We start our investigation of the phase structure of the lattice model with the case of zero 
temperature and zero chemical potential. To leading order in expansion, the ground 
state of the model is determined by the pair of saddle point equations 

^ cjL -6mL _ ^^^p^ + ^ 0, (3.1) 

= TTi f ^ - F(aL, t:l)\ = 0, (3.2) 



9t^l \gl 



where 



^^''^'''^^ (27r)2E.sin2e. + (^L + E.(l-cose.))2 + 7ri' ^^"^^ 

^ ^^ m 2E.(l-cose.) 

^KCJ^.^L) (27r)2E.sin2e. + (a^ + E.(l-cose.))2 + 7ri- ^^"^^ 

The second equation allows a non-trivial solution with ttl ^ 0, which corresponds to 

spontaneous breakdown of parity symmetry, in addition to the parity-symmetric solution 

TtL — 0. If the phase transition between the two phases is continuous, the phase boundary 

separating them can be determined by examining the limit of the parity-broken solution 

toward tt/, — > 0. This yields the pair of equations 

= <7lF(^l, 0) + G(ai, 0), (3.5) 
4 = ^(^L,0). (3.6) 



given 


by 










2 


d^VL 




1 




) - 




7TL=0 





The second equation ensures that the pion mass M,r vanishes along the phase boundary 

F{aL,0). (3.7) 

The structure of the phase boundary was previously studied for the case of an equal coupling 
gl = gl in Refs. We have to extend the analysis treating the couplings as independent. 

The pair of equations (|3.5| - |3.6|) defines a surface in the three-dimensional parameter space 
{g1, g'^^dniL). We shall tentatively call this surface as parity phase boundary. In Fig. ^ we 
show the section of the surface on the {6mL,gl) plane at g"^ = 0.8. The region above the 
curve corresponds to the parity-breaking solution ttl ^ 0, while vr^ = below it. As one 
moves along the curve from right to left, the value of ctl taken as the parameter of the 
curve decreases from positive to negative values. The bottom of the three valleys of the 
curve reaches g^ = at = ctlo = 0, —2, and —4 from right to left, respectively, where 
F{aL, 0) = +00 and the mass parameter 6mL equals 6mi = SmLQ^gl) = a^Q — g^G^aLo, 0). 

The horizontal line marks the position of equal coupling g"^ = g^. We observe that there 
are three parity-broken intervals in between the parity-symmetric ones. Toward the weak- 
coupling limit g^ ^ 0, the horizontal line moves toward the bottom of the valley, which in 
turn converges toward 6mLo{0) = aio = 0, —2, —4. Hence each of the parity-broken intervals 
narrows to a point at g^. = 0; the point at Smi = is the conventional continuum limit, 
while 6mi = —2 {6mL = —4) represents the point where the doublers with momentum 
^ = (0, vr) and (tt, 0) (^ = (vr,??)) become massless physical modes. These are the results 
previously obtained in Refs. PJ^. 

As we have emphasized in the previous section, however, a proper continuum limit has 
to be taken along the line specified by ( |2.27| - |2]28| ). The point T{g'^) corresponding to this 



line without 0{a) corrections, given in ( |2.33| ), is marked by a cross in Fig. |^. The parity- 
broken interval is significantly narrower near the point, and the location of the point itself 
is indistinguishable from the curve in the scale of the figure. A more detailed study of the 
region near the point is clearly needed. As we shall see, the phase structure near the point 
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is far more complex than would seem from Fig. |^. An indication is a self-crossing behavior 
of the curve already visible for the central valley in Fig. ^ which also occurs for the right 
and left valleys. 

In Fig. ^ we plot an expanded view of the region near the tuned point T{g'^)] we take 

= 2.0 for this figure as the structure is enhanced and so easier to draw. The thick curve 
represents the parity phase boundary, while the set of thin curves are lines of constant ttl, 
as determined from ( p.l[ - p.2|) , for < tt^ < 0.12 in steps of 0.02. A number of important 
features are revealed in this figure. First of all, the boundary curve crosses with itself at the 
point B (this point differs from the tuned point T((7^) marked by a cross, to be discussed 
below). To understand why the crossing takes place, we note that the right-hand side of 
( p.5|) behaves as (0"/, ln(T|)/27r + 2Ci near ctl = 0. Hence ( p.5|) regarded as an equation of 
aL for a given 6171^ has three solutions in the neighborhood of 6171^ = —2Cig'^, leading to 
the self-crossing behavior. Since the existence of three solutions is valid for arbitrarily small 
values of g^, the crossing persists down to gl = 0. Quite clearly there is no guarantee that 
the boundary curve represents the true phase boundary below the crossing point. 

Another important feature is that the lines of constant ttl spill beyond the boundary 
curve to the right of the crossing point. In the triangular region DBE formed by the two 
segments of the boundary curve DB-BE and the envelope ED of lines of constant ttl, the 
saddle point equations ( p.l[ - |3.2|) have three solutions, one with tt^ = and two more with 
7Tl ^ 0. Which of the solutions represent the true minimum can only be determined by 
examining the value of the effective potential. This is a typical situation where a first-order 
phase transition occurs. 

A detailed numerical investigation of the effective potential shows that the true phase 
boundary runs as follows. Starting from right, the true phase boundary coincides with the 
parity phase boundary, and hence being of second order, from the point A to the point D 
where the envelope of lines of constant vr^ starts out. Beyond this point, the phase transition 
becomes of first order, and runs along the dashed curve from the point D to the point F 
which is located on the parity phase boundary curve below the crossing point B. At the point 
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F, the phase transition returns to second order, and runs along the boundary curve from 
the point F through the crossing point B toward the point C. Altogether the parity-broken 
phase occupies the region above the line ADFBC. 

Below this line the ground state of the model can be examined by the effective potential 
with 7iL = 0. The potential no longer depends on g"^, and as a function of ai has a double 
well structure. As one varies 6mL below the point F from right to left, the true minimum 
switches from a positive to a negative value of ai discontinuously. This leads to a first-order 
phase transition along the line FH reaching down to = 0. We shall call this line as the a 
phase boundary. Qualitatively, one may describe the behavior of the system near and below 
the point F as "almost" exhibiting spontaneous breakdown of chiral symmetry. 

We emphasize that the intricate structure involving first-order transitions described 
above originates from the 0{a^) term of the effective potential Vl {i.e., 0{a) in physi- 
cal units) which contains terms cubic in ai- This contrasts with the leading O(a^) terms 
in ( |2.26| ) which gives rise to only simple second-order transitions on the {SrriL, g^) plane. 
Indeed, including the 0{a^) contributions given by 

5Vl = —Csal + 2 k - ^ ln(ai + Tri)! a^ai + Trf) (3.8) 
o L ovr 

with 

C3 = 0.0647- ■■, Ci = 0.0904 ••• , (3.9) 
to the effective potential ( p.26|) , one quantitatively reproduces the features of the phase 



diagram shown in Fig. 0. 

The phase structure on the {6mL,gl) plane we have discussed for g1 = 2 remains the 
same toward (7^ — except that the structure as a whole moves down toward g^ = Q 
and becomes narrower horizontally. In the weak-coupling limit g^. = 0, the boundary line 
degenerates tol/gl = F{6mL, 0) which does not have a self-crossing. The first-order line DF 
disappears, and the entire phase boundary becomes of second order. We thus find the phase 
structure in the three-dimensional parameter space {g"^, g'^, 61711) as schematically drawn in 
Fig. |. 
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Let us now discuss the continuum limit, in particular, how the phase structure relates 
with the tuning of parameters specified by ( |2.27| - |2?28D . The tuned point Tig"^) without 0{a) 
terms (cross in Fig. |^) is located slightly to the left and below the point B in the parity- 
symmetric phase. When one takes the coupling g"^ to zero, the relative location of the tuned 
point and the a phase boundary remains the same. Therefore taking the continuum limit 
exactly along the tuned point Tig"^) leads to a negative value of cr. This, however, does not 
mean that a continuum limit with a positive value of a is not possible. 

As we have already pointed out, the structure around the point F stems from the (9(a^) 
terms of the effective potential, and hence its size shrinks as a power of a. One can estimate 
the effect of adding (|3^ ) to ( |2.26| ) by an order counting. Since 6171^ is the coefficient of the 
linear term in ( p.26| ), the 0{a^) correction (|3.8|) causes a shift of form 

Stul ^ SrriL + 0{al) ^ 6mL + 0{a^) . (3.10) 

Hence the structure has a size of 0{a'^) in Sttll. By a similar argument for g^, we conclude a 
size of 0{a) for g^. The point of tuning, on the other hand, has a degree of freedom of 0{a'^) 
for 6m L and 0{a) for g"^ . Therefore one may shift the point of tuning within a rectangular 
region of size O(a^) times 0{a) around the point T((yf^), including a shift to the right of the 
point F, in which case the continuum value of a becomes positive. 

In order to demonstrate that the correct continuum limit is obtained within the freedom 
of the choice of the parameters described above, we examine the scaling of two physical 
observables a/A and (M^r/A)^ along the lines given by 

SrriL = -2Cigl + Am(aA)2, (3.11) 
^ = 4 + + A(7^aA , (3.12) 

where Am and Ag^^ are 0{1) constants. We show the solution of the saddle point equation 
( p.l| ) for cr/A = o"i/(aA) for Am = 0, 1,2 and Ag^^ = in Fig. ^ as a function of aA = 
cexp{—7r/gl) where A is defined in ( p.29| ). One sees that a/A converges to the correct 
continuum result a/A = 1 whether the limit is taken on the left (Am = 0) or the right 
(Am = 1,2) of the point F of the phase boundary. 
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The pion mass is more interesting since the mechanism [^,0 which ensures = 
even for finite a no longer works along the segment DF in Fig. ^ because of the first-order 
nature of the phase transition in this part of the phase boundary. We select Am = 1, 1.5, 2.0 
and tune A^f^r so that the point (5mL, g^) given by ( |3.11[ - |3.12D is placed just under the first- 
order line DF in the parity-symmetric phase. Solving the saddle point equations ( ^.l[ - p3.2D for 
ai and substituting the result into (|3J| ), we obtain the curve for (Mt^/A)^ plotted in Fig. ^ 
as a function of aA. The squared pion mass vanishes as a power of a toward the continuum 
limit, demonstrating that the first-order nature of the parity-breaking phase transition along 
DF does not cause problems. 

So far we have concentrated on the phase structure near Sttil = 0. The phase structure 
around 6m ^ = —4 is obtained by a refiection 6m ^ —4 — 6mL. The phase structure near 
6mL = —2 is different. It is symmetric across the 6mL = —2 plane, and the parity phase 
boundary turned out to be always of second order. In addition a a phase boundary of first- 
order transition across which the sign of {a^ + 2) fiips exists at 6mL = —2, similar to the 
line FH in Fig. |. 

In the following we concentrate on the phase structure around 6mL = relevant for the 
usual continuum limit. 



IV. PHASE DIAGRAM AT T / AND fi = 

Our analysis of the phase diagram for the case of finite temperature and zero chemical 
potential proceeds in essentially the same way as for the zero-temperature case; choosing 
a temporal lattice size Nt and a value of g"^, we analyze the solution of the saddle point 
equations and the value of the effective potential given by ( |2.11[ - |2rT5|) with fii = as a 
function of g"^ and 6mL. 

In Fig. 1^ we show the result of such an analysis at g"^ = 2 for a set of temporal lattice 
size Nt- For large Nt, the topological structure of the phase boundary remains the same 
as at Nt = oo corresponding to zero temperature. In particular the left-half of the phase 
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boundary hardly changes. On the other hand, the right-half of the boundary with the first- 
order segment moves toward upper left. The lower end point of the first-order segment slides 
up along the second-order part of the boundary for Nt = oo- The vertical a phase boundary 
also moves toward left with the end point. 

This structure changes at a threshold value of N^, which equals = 13 for g"^ = 2 
chosen for Fig. ^ when the first-order segment of the boundary disappears. The entire phase 
boundary becomes a smooth line of second-order transition, along which ai continuously 
varies from positive to negative values as one traverses the boundary curve from right to 
left. As may be inferred from this, the vertical a phase boundary also disappears at the 
threshold temporal lattice size. 

To the extent that the first-order a boundary is interpreted as indicative of spontaneous 
breakdown of chiral symmetry, its disappearance may be regarded as signifying restoration 
of chiral symmetry. Indeed the value of ai becomes small when the a phase boundary 
disappears. We illustrate this point in Fig. ^ which shows how the effective potential changes 
from a double-well to a single-well structure across the threshold value of Nt at the tuned 
point Tigl) for gl = 0.8. 

The structural change found as a function of Nt for a fixed g^ can be translated into 
that for a varying g"^ with a fixed temporal size Nt '■ as one lowers g"^, there is a threshold 
value g'^ = g'^^.i^T) below which the first-order segment of the parity-broken phase boundary 
and the vertical first-order a boundary cease to exist. This leads to the three-dimensional 
phase structure for a fixed and finite Nt schematically drawn in Fig. p. Contrary to the 
Nt = oo case, the boundary curve at g"^ = does not reach the value g"^ = since it is given 
by l/gl = F{6mL,0) where the function F{SmL,0) computed for a finite Nt is finite and 
positive at its maximum at 6171^ = 0. 

Let us consider the continuum limit of the model at finite temperature. This limit is 
defined by simultaneously taking A''^ 00 and gl such that 



T 
A 



exp(7r/^g) 
cNt 



(4.1) 
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is kept fixed. The coupling g1 and the mass parameter SrriL have to be tuned according to 
( p.27| - |2^2^ ) to restore chiral symmetry. An alternative way to take the limit, which is closer 



to the method practiced in numerical simulations, is to vary g1 for a fixed Nt- For each 
this yields physical observables as a function of T/A, and the continuum limit is obtained as 
a limit of this function as Nj- oo. We follow the latter approach in the following analysis. 

A qualitative picture of how finite-temperature transition appears in this approach is as 
follows. The continuum limit is taken along the line satisfying (p.27| - p^28D , which converges 



to {gl, gl,6mL) = (0,0,0). Above the threshold value gl = g^^^Nx), this line runs along 
either side of the lower edge of the parity-broken phase where physical quantities take values 
similar to those at zero temperature. In particular ai has a non-zero value of positive or 
negative sign depending on whether the line is on the right or left of the vertical a phase 
boundary. However, as g^ is lowered below the threshold value, the parity-breaking phase 
boundary moves up toward larger g"^, and the a phase boundary disappears. Once the 
line enters into this region, the value of cr^ becomes small. In other words, for each fixed 
temporal size, the threshold value of g"^ marks the region of finite-temperature transition. 

We illustrate the description above in Fig. |TD| where we plot cr/A as a function of T/A 
calculated along the line T{gl) without 0{a) corrections for several values of Nt, together 
with the curve in the continuum limit. We observe that the curves for finite Nt smoothly 
approaches the continuum result drawn by a thick line as Nt increases toward oo. 

Several remarks are in order with this figure: (i) The negative sign of a is due to the 



fact, already remarked in Sec. p| , that the continuum limit is taken on the left side of the 
a phase boundary, (ii) In the region close to the continuum critical point T/A ^ T^/A = 
exp(7£;/7r) = 0.56694- ■ -, o"(T) has a small discontinuity, as shown in the inset of Fig. [1^, 
which reduce with increasing Nt- This is due to the fact that the tuned point T((7^) goes 
through the vertical a phase boundary as g^ is decreased just before the threshold value is 
reached (see the relative location of the boundary of Nt = 13 and T{g'^) in Fig. |^). (iii) The 
value of cr(T) after the discontinuity is very small but non-zero for a finite Nt- It exactly 
vanishes for T >Tc only in the limit of Nt oo. 
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Another possible choice for the continuum hmit is to run along the bottom of the valley of 
the parity-breaking phase boundary for each value of A^^ (line EJB in Fig. The limiting 
point B of the line at = has a non-zero value of g^. This point, however, reaches 
the correct continuum limit {g^, gl,6mL) = (0,0,0) as Nt oo. Above the threshold 
value of gl (point J) where the a phase boundary exist, <j(T) has two solutions of opposite 
sign(a"^, a^^). At the bottom of the valley, these two solution are the balancing minima of 
the effective potential, i.e., 

Vdat) = VLial) . (4.2) 

At the threshold value of g"^, cr^ and aj^ merge through a second-order singularity corre- 
sponding to the second-order phase transition (line JI) marking the end of the first-order a 
phase boundary. 

With this choice of the continuum limit, we obtain the behavior of o""'"(T) and cr~{T) 



shown by dashed lines in Fig. |Tl] and Fig. |T2|, exhibiting convergence to the continuum result 
(thick line). A notable property is that the system exhibits a second-order phase transition 
already for finite values of Nt as the line passes through the point J in Fig. |^. Also shown 
in Fig. |TT| by dotted lines is the difference of the two solutions divided by 2, {a^^ — al)/2, 
which shows a faster convergence to the continuum result. 

V. PHASE DIAGRAM FOR / 

A. The Case of /x / and T = 

We now consider the phase diagram for non-zero chemical potential, starting with the 
case of zero temperature {Nt = oo). In this case, from (|2.11[ - ^rT5D , the effective potential is 
given by 

Vl = TT^i^^L - SniLy + —ttI - / \n{aL + 2 - cos^ + max(E, /i^) . (5.1) 
Zg^ Zg^ j-ttZti 
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As one expects from this expression, as long as /j^l is small, the phase boundary on the 
{6mL,g'^) plane does not move from that with fi^ = 0. However, as we show in Fig. |13 



for gl = 2, when fii exceeds a threshold value {fiLc ~ 0.09 for Fig. |T^, the vertical a 
phase boundary splits into two vertical first-order boundaries such that ai between the 
two boundaries is near zero. The end points of the two boundaries slide upwards in the 
opposite directions along the parity phase boundary for Nt = oo, and the phase boundary 
connecting the two end points, which is of first-order, is convex and moves upward toward 
larger g^. This behavior reflects the fact that the effective potential, which has a double 
well structure at /i/, ~ develops an additional minimum close to o"l = as /i^ increases 



beyond the threshold value. See Fig. Which of the three minima represents the true 
minimum depends on the value of 6m l, which leads to the double first-order boundaries 
seen in Fig. |13[ 

An important point with the behavior described above is that the two vertical first-order 
boundaries move away from the edge point F of the parity phase boundary at = 0. Since 
the point ( |2.27| - P^^5D for taking the continuum limit is located near the point F, this means 
that the system, initially in one of the phases with a non-zero value of o"/, at /i/, = undergoes 
a first-order phase transition when the boundary moves over the point as fi^ increases. After 
this takes place, the system is in the middle phase with ctl ~ 0, which qualitatively means 
that chiral symmetry becomes restored. This, then, is the mechanism by which the lattice 
model yields the first-order finite density phase transition at zero temperature. These points 
are illustrated in Fig. |1^ for g"^ = 2.0. For the tuned point T((yf^) , for example, the left vertical 
phase boundary moves over the point between fii = 0.12 and 0.13. In Fig. |16| we translate 
the behavior into a schematic three-dimensional phase diagram for A^^ = oo and a fixed 
value of fiL. 
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B. The Case of /x 7^ and T ^0 



Finally we consider the finite-density transition at finite temperature. We examine this 
case by decreasing the temporal lattice size from = 00. The phase structure for A^^ = 00 
found in Sec. |V A| implies that the mechanism of generation of a first-order transition through 
a splitting of the a phase boundary as fii increases will remain operative for large temporal 
size Nt- We found this to be the case as shown in Fig. 17 for Nt = 120 at g"^ = 1.5. 

On the other hand, for small temporal sizes, we expect the behavior to become similar 
to the finite-temperature and zero chemical potential case examined in Sec. |TI|. We plot 



an example of this case in Fig. |I8| for which Nt is decreased to 28 : the parity-broken 
phase shrinks, and the first-order segment of the phase boundary and the a phase boundary 
disappear as one increases fiL. In this case the phase transition is of second order in the 
continuum limit. 

We plot in Fig. |19| and Fig. ^ curves of (j(T, /i)/A calculated as a function of /i/A for a 
fixed value of T/A along the tuned point T((7^). A clear discontinuity observed in Fig. ^ 
for T/A = 0.2 contrasts sharply with the convergence to second-order behavior seen for 
T/A = 0.5 in Fig. This difference reflects the change of order of flnite-density transition 
for large and small temporal lattice sizes discussed above. In the continuum limit we thus 
flnd a line of phase transition on the (/i, T) plane which switches from being of flrst-order 
to second order as the critical temperature increases. This agrees with the phase diagram 
of the continuum theory drawn in Fig. |l|. 



VI. CONCLUSIONS 

In this article we have investigated the phase structure and the continuum extrapolation 
of the two-dimensional lattice Gross-Neveu model at flnite temperature and chemical poten- 
tial. The choice of the Wilson action for fermions leads to the existence of a parity-broken 
phase separated from the parity-symmetric phase by a two-dimensional boundary surface 
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in the three-dimensional space of the couphngs {g^, gl,SmL) that specifies the modeL The 
phase transition across the surface is generally of second order. However, in the region rel- 
evant for the continuum limit, 0{a) effects of the Wilson term gives rise to a complicated 
structure, rendering the transition to be of first order in a region whose size vanishes in the 
continuum limit. The parity-symmetric phase is also divided into subregions by a first-order 
phase boundary corresponding to different minima of the effective potential. We found that 
these detailed structures intertwine to lead to the continuum phase diagram reported in 
Ref. [|11|. 

An interesting issue is how much of the detailed structure we found for the Gross-Neveu 
model has relevance for QCD in four dimensions. Of particular interest is the question that a 
part of the parity-flavor breaking phase boundary may be of flrst order. While the necessity 
of two independent couplings and gl appears to be a speciflc feature of the Gross-Neveu 
model, arising from the presence of four-fermion couplings in the model, a distortion of the 
effective potential for effective meson flelds in the a direction by terms odd in a would be 
present also for QCD, possibly triggering a flrst-order phase transition. 

The relevance of the structure we found for flnite chemical potential is more likely; the 
emergence of three regions corresponding to three different values of a for large /i is a 
physical effect, not speciflcally tied with the use of the Wilson actions for fermions. The 
crossover from a flrst-order to second-order transition depending on the critical temperature 
means, however, that the order of the flnite density transition may depend sensitively on 
the parameters of the model, perhaps the space-time dimension playing an important role 
rT|. These, then, are the issues which have to be directly addressed in QCD. 
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FIGURES 
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FIG. 1. Phase diagram of the continuum Gross-Nevcu model in the (///A, T/A) parameter 
space. SoUd and dashed curves represent a second- and a first-order segment of phase transition 
restoring chiral symmetry. 
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FIG. 3. Magnified view of the region near the tuned point T{gl) for gl = 2.0. Thick 
curve represents the parity phase boundary while thin fines are contour curves of fixed ttl val- 
ues (ttl = 0.02, 0.04, 0.06, 0.08, 0.10, 0.12). True phase boundary consists of AD (second order), 
DF(first order) and FBC(second order). 
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FIG. 4. Schematic phase diagram for Nt = oo and fXL = 0. The phase bound- 
ary of the parity-broken phase AOB-CDEF forms a blade whose edge touches the point O 
(S'ctjS'^) '^"^i) = (0)0)0)- The first order part EOD shrinks to a point at g'^ = 0. Parity-symmetric 
phase is divided into two phases with ctl < and ctl > by a first-order a phase boundary DOG 
which vertically drops from the edge OD. 
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FIG. 5. Order parameter a in units of A as a function of lattice spacing aA = cexp(— Tr/g"^) 
calculated for three choices of continuum extrapolation specified by Am for = 0. Sign of a is 
negative for Am = 0. 
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FIG. 6. Pion mass squared as a function of lattice spacing aA under continuum extrapolation. 
See text for tuning of couplings employed. 
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FIG. 7. Temporal lattice size dependence of phase structure plotted for Nt = oo, 16, 14, 13, 
12, 10, 8. Phase transition is of first order along dashed lines, while it is of second order along solid 
lines. 
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FIG. 8. Effective potential for = 0.8 for temporal lattice size Nt = 160 and Nt = 154. 
is tuned according to T{g^). 
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FIG. 9. Schematic phase diagram for fixed and finite Nt- For smaU g'^, first-order par- 
ity-breaking bmHidary (FJE) and a phase boundary (EJIH) disappear. 
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FIG. 10. Order parameter a as a function of temperature T for = evaluated for several 
temporal lattice size along the point T{g'^) without 0{a) corrections. Inset shows an expanded 
view of the critical region. 
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FIG. 11. Positive solution as a function of T for fx = 0. Couplings are tuned along the edge 
of the parity-breaking phase boundary (line EJB in Fig. ^) for each value of Nt ■ The combination 
(cj^ — cr^)/2 is also plotted (see curves specified as "(mid.)"). Inset shows an expanded view of 
the critical region. 
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FIG. 13. Phase structure for = 0.110 and Nt = oo at = 2.0. Hatched region is the new 
phase with cr^ 0. Phase transitions between this phase {ttl = 0, (Ti 0) and other three phases 
{ttl 7^ 0, (Tl < 0, aL > 0) are all first order. 
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FIG. 15. Same as Fig. |T| for m = 0, 0.10, 0.11, 0.12, 0.13. Symbol T represent the tuned 
point T{gl). 
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FIG. 16. Sketch of phase structure for fixed and finite jiL for Nt = oo. Below a threshold 

value of the vertical a phase boundary splits into two, both being of first order, enclosing a 
region with aL ^ 0. 
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FIG. 17. Effect of /il on the pliase structure for large A^r- 
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FIG. 18. Effect of /x^ on the phase structure for small Nt- 
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FIG. 19. —a as a function of for T/A = 0.2. Couplings are tuned at T{g^). 
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FIG. 20. Same as Fig. |T| for T/A = 0.5. 
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